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پایش سطح زیرکشت برنج با استفاده از سیستم اطلاعات حغرافیایی و تکنیک سنحش از دور 


۲ ۱ ۲ ۱۳ 
ارش مصری . فاطمه رحیمی احدادی 5 ایرج باقری 


تاریخ دریافت: ۱۳۹۹/۱۲/۱۰ 


تاریخ پذیرش: ۱۴۰۰/۰۴/۰۹ 


چکیده 


سنجش از دور یکی از ابزارهای کارآمد برای بررسی روند تغییرات سطح زیر کشت محصولات کشاورزی و باغی در سطوح وسیع و زمان کوتاه است. 
منیاست‌گذاران با اگاهی از این اطلاعات: م تو اند کضیمات. صحیم و یموق قاشقه باشتت مظالدی حاضن با هدف کین بطم زیر کشت شالیزازهاق 


برنج در بخش کیاشهر استان گیلان انجام شد از تصاویر سنچنده 1111 ماهواره لندست ۵ و سنجنده .01 ماهواره لندست ۸ به‌منظور تهیه نقشه‌های 
کاربری اراضی استفاده شد. ابتدا؛ تصحیح هندسی و اتمسفری بر روی تصاویر صورت گرفت. سپس با استفاده از الگوریتم طبقه‌بندی نظارت شده حداکثر 
تیال ها اریز ازاخی مه با هفت کاربری هام اراس ی نگل ساره خن تک سای وان ممعاطی ای »نتفای 
ماسه‌ای و سایر اراضی تهیه شد. در ادامه. مساحت هر یک از کاربری‌ها محاسبه شد و روند تغییرات» مورد مقایسه قرار گرفت. دقت کلی و ضریب کاپای 
طبقه‌بندی به‌ترتیب معادل 8۸/۴۵ و ۰/۹۸ برای سال ۸۲۰۰۰ ٩۷/۵۹‏ و ۰/۹۷ برای سال ۲۰۱۰ و ۸/۷۲ و ۰/۹۸ برای سال ۲۰۲۰ به‌دست آمد. 
نتایج نشان داد که اراضی برنج در یک بازه ۲۰ ساله با کاهش ۶/۹۴ درصدی همراه بوده به‌طوری که مساحت آن از ۱۱۰۸۰/۶۶ هکتار در سال ۲۰۰۰ به 
۹ عهکتار در سال ۲۰۲۰ رسیده است. همچنین» در این مدّت مناطق مسکونی و جنگل‌های تنک به‌میزان ۶۷/۹۴ و ۱۸/۷۳ درصد رشد کرده‌اند. 
اما جنگل‌های نیمه‌انبوه» مناطق آبی و پهنه‌های ماسه‌ای به‌ترتیب ۱/۳۲ ۴/۹۱ و ۶۱/۴۸ درصد کاهش داشتند. با توجه به نتایج توجه جدّی به تغییر 


کاربری اراضی برنج و تخریب جنگل‌ها ضروری می‌باشد. 


واژه‌های کلیدی: حداکثر احتمال جنگل تنک طبقه‌بندی نظارت‌شده. کاربری اراضی 


مقدمه 

برنج چهارمین غله پرمصرف جهان و غذای بیشتر از نصف مردم 
به‌خصوص در کشورهای در حال توسعه است ,۲(۵۲۷15۳072060) 
1010101 :2019 ,۳۲۸۵ :2011 ,۳۲۱200871 عک مصه۷]21 
(2017 ,0۵1 ک ,۹2125 ,ورتتاطل ۰0۷ همچنین اين محصول یک 
سوم از کل سطح زیر کشت جهانی غلات را در برداشته که حدود ۳۵ 
( ع ,920727طوظ تفه ردهمنام عصفوظ . تمصم‌مهکز 
7 ,25111). بیشتر از ٩۰‏ درصد محصول برنج جهان به‌دلیل آب 
و هوای گرم و مرطوب (عرض جغرافیایی ۵۵ درجه شمالی تا ۲۶ 
درجه جنوبی) در قاره آسیا تولید می‌شود ( ,نصطی11هط06۵۱ بطانز۸ 
7 ,۵8866۲561۷25 عک رصع61 ۳2۳۵۵۵1۷ ,ط24اع162). سطح 


که بیشتر این اراضی در استان‌های شمالی» شامل استان‌های مازندران 


۱- دانشجوی کارشناسی ارشد گروه مهندسی بیوسیستم دانشگاه گیلان» رشت» 

ای 

۲- استادیار گروه مهندسی بیوسیستم» دانشگاه گیلان. رشت ایران 

(*+- نویسنده مسئول: ۴ ۵ 2 _تصنطه: :اتفصط) 
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با ۲۸/۵ درصد. گیلان با ۳۱/۹ درصد و گلستان با ۹/٩‏ درصد قرار 
دارد. همچنین خوزستان با ۹/۲ درصد و فارس با ۲/۸ درصد در 
رده‌های بعدی قرار دارند > بلعذ2 ,طع220ن2ک ,۲2200056) 
(2013 ,نطلاه0ط۸. زندگی روستاییان و معیشت آن‌هاء به‌ویژه در 
نواحی شمالی کشور به کشاورزی وابسته است. استفاده بهینه از 
زمین‌های زراعی» بذرهای اصلاح شده ماشین‌آلات» سیلوهاء انبارها و 
غیره نیازمند مدیریت صحیح و دقیق مسئولان کشوری و استانی 
است. سیاست‌گذاران» با در اختیار داشتن سطح زیرکشت محصولات 
می‌توانند خدمات و تحهیزات را منطبق با داشته‌های موجود ارائه کنند. 
همچنین اين اطلاعات پتانسیل واقعی نواحی کشاورزی را مشخص 
کرده و از تغییرات کاربری آن جلوگیری خواهد کرد ( صلعه2 
,۱09060 ماو ی رتلممطنظ ههد رن۲۱۳0۱۵2۵0۵ 
89 با توجه به روند رو به رشد جمعیت و افزايش تغییرات کاربری 
اراضی در بخش کشاورزی نیاز به اطلاعات دقیق و به‌روز از سطح 
زیرکشت محصولات. میزان عملکرد محصول در واحد سطح و بررسی 
تغییرات آن‌ها جهت مدیریت صحیح منابع اجتناب‌نایذیر است. 
به‌دست‌آوردن اين اطلاعات با استفاده از روش‌های سنتی» زمان 
زیادی می‌برد و صرفه اقتصادی ندارد > ,12:6 ,طعصن۹ ,۳۳۵۵۵۵) 
(2009 ,.۵ ۶۶ حعطگ :2007 ,وماع]ه>1. در ایران» لا تعیین 


مصری و همکاران. پایش سطح زیر کشت برنج با استفاده از سیستم اطلاعات جغرافیایی ... ۵۱۷ 


سطح زیرکشت محصولات با سه روش کارشناسی» فهرست‌برداری و 
استفاده از فناوری‌های جدید انجام می‌گیرد. تخمین سطح زیر کشت در 
روش کارشناسی دارای دقت بسیار اندکی است و نمی‌توان از نتایج آن 
در تصمیم‌گیری‌های مهم کشور استفاده کرد. در روش برآورد از طریق 
فهرست‌برداری» هرچند سطح زبرکشت محصول تخمین زده می‌شود. 
نیروی انسانی بسیار بالا خواهد بود. سومین روش» استفاده از 
فناوری‌های نوینی چون سنجش از دور و سامانه اطلاعات جغرافیایی 
است (2009 ,.۵ ۶ ۳۲02۵0201 صمتع212). این روش دارای 
محاسن فراوان بوده که در روش‌های سنتی قبلی وجود ندارد. علاوه 
یعنی هزینه و زمان کمتر تهیه اطلاعات را نیز برآورده می‌سازد 
,12600 ,وعله ۷۷ :2016 مصقصصطهای عک روتاهاصتا۱۷]۵ رتحصع/قطک) 
(2020 ,10۷611167 ک. در دهه‌های گذشته استفاده از سنحش از دور 
در بخش کشاورزی افزایش یافته و اکنون سنجش از دور به ابزاری 
موثر 9 پرکاربرد برای نقشه برداری از محصولات تبدیل شده است 
بااع050ظ مطم۱2۷۱۵۵2۵0 مصفطیف :2013 رته‌عتهها2اض) 
(2019 ,۱6۱908 »6 ,۳606. تعیین عملکرد محصولات. تمایز میان 
هستند »4 ,2206 ,1۳08۵ مقطعصنه :2014 ,اه ۶ ۷۵8021) 
(2019 ,2120. در ادامه به برخی از نتتایج مطالعات انجام شده در 
زمینه استفاده از تصاویر ماهواره‌ای برای تهیه نقشه‌های کاربری 
اراضی و تعیین سطح زیر کشت محصولات کشاورزی پرداخته شده 
است. 

در پژوهشی از تصاویر سنجنده ۸۷۳۲ ماهواره 0۵۸۸۵ 
به‌منظور طبقه‌بندی محصول برنج و برآورد سطح زیرکشت آن در 
استان گیلان استفاده شد. برای این منظور از روش طبقه‌بندی چند 
زمانه متکی بر استخراج روزانه ۱۲۲۷1 بيشینه به منظور حذف ابر 
استفاد شد. نتایج نشان از دفت ٩۱/۹۶‏ درصدی سطح زیرکشت 
محاسبه شده در مقایسه با نقشه مرجع دارد ( ک تامص۸ تتحوحم 
1 ,۹2720 4مصصعطم‌صنآض). در پژوهشی دیگر برای تعیین 
سطح زیرکشت دو محصول برنج و سویا در چهار حوزه آبخیز استان 
گلستان از سه روش طبقه‌بندی نظارت شده استفاده شد. نتایج نشان 
داد که روش‌های طبقه‌بندی حداکثر احتمال" و کم‌ترین فاصله از 

بر رنه و : ۲ 

میانگین به‌ترتیب برای تخمین سطح زبرکشت برنج و سویا روش 
مناسبی هستند (2019 ,تحطع12 ی رتقاجصهک رناذ۷ ۱۷۵ تاطافدن/ 
نتایج تحقیق دیگری حاکی از آن بود که سطح زیر کشت محاسبه 
شده با استفاده از دو روش طبقه‌بندی حداکثر احتمال و شبکه عصبی 
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به‌ترتیب با اختلاف ۱۶/۸ و ۱۳/۲ درصدی می‌تواند مساحت 
محصولات مختلف را تخمین بزند ( بتعحط‌اطهم ,تناممتلش 
4 ,۹006 عک بلت۲2-تا0مدعدداش). در پژوهشی دیگر به 
منظور برآورد سطح زیرکشت باغات از داده‌های 1۳5-۳6 استفاده شد. 
طبقه‌بندی داده‌ها با الگوریتم‌های مختلف و به دو صورت هفت 
کلاسه (طبقات کاربری اراضی) و دو کلاسه (باغ و غیرباغ)انجام شد. 
نتایج نشان داد که الگوریتم حداکثر احتمال در طبقه‌بندی به صورت 
هفت کلاسه و دو کلاسه دارای بیشترین صحت کل و ضریب کاپا" 
بودند (2019 ,تهتقطگ ع رنصهک نهمحصطم ,تععصناه ۷). در تحقیقی 
رابطه رگرسیونی بین چندین شاخص پوشش گیاهی (۸۱00۷1 
۲ ۷۲ و ۷۲ استخراج شده از تصاویر ماهواره لندست ۷ 
و عملکرد برنج در سال‌های ۲۰۱۲ و ۲۰۱۳ ارزیابی شد. مقایسه 
مدل‌های تولیدی نان از برفری شاخص ۲۷۲ داشت 
(2018 ,تلدطمطنوه]۱۷ عک ,تن بهتنعدط ,تانامطع۷2). تصاویر 
چند زمانه ماهواره لندست *1 به منظور تعیین سطح زیرکشت و 
سن گیاه برنج در منطقه بالی اندونزی استفاده شد. برای این منظور 
شاخص جدیدی به نام 0۷1 توسعه داده شد ( ,منطعئل ,دکته۵ 
0 ,۳۲0920 4). در پژوهش دیگری نتیحه شد که از بین 
شاخص‌های پوشش گیاهی, شاخص 1360۷1 بیشترین همبستگی را با 
عملکرد برنج داشته و با آن می‌توان عملکرد شالیزارهای برنج طی ۲۰ 
تا ۶۰روز قبل از برداشت را پیش‌بینی کرد ( ع مقطهگ بلنا 
0 ,۳]6۳00۲2۳). تحقیقی به منظور شناسایی محصولات با 
استفاده از تصاویر ماهواره 5 ٩۳0۲‏ انجام گرفت. در این تحقیق از 
پنج روش طبقه‌بندی برای شناسایی انواع محصولات استفاده شد. 
نتایج نشان داد که دو روش طبقه‌بندی حداکثر احتمال و ماشین‌های 
پشتیبان بردار " از دیگر روش‌های مطالعه شده بهتر هستند ( ,۷۵۵۵ 
1 ,۱۷۲0۲06 ک باانته۳۷). در تحقیقی دیگر از تصاویر ]7۷ 
ماهواره لندست به منظور طبقه‌بندی کاربری اراضی استفاده شد. در 
این تحقیق از روش‌های نظارت‌نشده و سیزده مورد از روش‌های 
نظارت‌شده برای طبقه‌بندی استفاده شد. نتایج نشان داد زمانی که 
نمونه آموزشی به تعداد مناسب انتخاب شود اکثر الگوریتم‌ها عملکرد 
خوبی دارند اما زمانی که تعداد نمونه‌های آموزشی کاهش می‌یابد. 
اختلاف دقت میان الگوریتم‌ها افزایش می‌یابد ,۷۷20 ,۷۷2082 ,1) 
(2014 ,ع005 عک ,تا در پژوهشی از داده‌های چند زمانه سنجنده 
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آآ) ماهواره لندست ۸ برای تعیین بهترین شاخص استخراج آب به 
منظور نظارت بر مناطق کشت برنج در ژاپن استفاده شد. پس از 
مقایسه شش شاخص, شاخص ۲21۷۲-051" که از باندهای سبز و 
مادون قرمز با طول موج کوتاه استفاده می‌کند بهترین عملکرد و 
کمترین میزان خطا در برآورد مناطق کشت برنج را داشت 
(2018 بلنایتطع1 عک ,مامحصقام بعناعد۳ع٩‏ رمامصیاهم) 

در سال‌های اخیر, پايش وضعیت سطح زیر کشت برنج. به‌عنوان 
یک محصول استراتژیک در استان گیلان» اهمیت بیشتری پیدا کرده 
است. دلیل اين امر, رشد مناطق ساحلی و توریستی استان و مهاجرت 
بی‌رویه و کنترل نشده ساکنین استان‌های جنوبی کشور به اين مناطق 
می‌باشد که اخیرا به‌دلیل پدیده خشکسالی در جنوب رشد بیشتری 
داشته‌است. این قبیل مهاجرت‌ها به شدت سبب تغییر و تبدیل اراضی 
کشاورزی به مناطق مسکونی شده است. علاوه بر اين» شیوع پدیده‌ی 
ویلاسازی» چه توسط شهرنشین‌های بومی استان و چه توسط اهالی 
استان‌های پرجمعیت و صنعتی اطراف منجر به تغییر کاربری بسیاری 
از اراضی» به‌خصوص اراضی شالیزاری توسط روستاییان شده است. 
روند این تغییرات در سال‌های اخیر به حدی چشمگیر بوده است کد 
کیان به یکی از یاو اصلی اشتاهای الکو شاه انگ: 
به‌طوری که توجه بسیاری از مسئولین استان را به خود معطوف کرده 
است. واضح است» قدم اول در مدیریت و کنترل معضل به‌وجود آمده 
امکان رصد وضعیت موجود و بررسی روند این تغییرات در چند دهه‌ی 
اخیر می‌باشد. هدف اصلی این مطالعه. بررسی سطح زیر کشت برنج 
قرع یک فاسله تما سایق سا ها ۵۱۳ ۱۳۹۹ ور 
بخش کیاشهر از شهرستان آستانه اشرفیه در استان گیلان با استفاده 
از فناوری سنجش از دور و سیستم اطلاعات جغرافیایی و مقایسه این 
تغییرات است. گفتنی است. این بخش یکی از مناطق ساحلی و 
گردشگری استان گیلان بوده و به جهت وجود پارک ملی بوجاق 
به‌عنوان یکی از زیستگاه‌های ارزشمند حیات وحش حائز اهمیست 
بین‌المللی می‌باشد. از اهداف فرعی تحقیق بررسی تغییر کاربری سایر 
اراضی» شامل اراضی جنگلی» مسکونی» آبی و پهنه‌های ماسه‌ای و 
میزان تبدیل آن‌ها به یکدیگر می‌باشد. 


مواد و روش‌ها 


بندر کیاشهره شهری در فاصله ۱۷ کیلومتری آستانهاشرفیه است 


در طول جغرافیایی ۲۲/۶۴ ۵۳ ۳۹ تا ۱۵/۲۶ ۱۱ ۵۰ شرقی و 
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عرض جنرافیایی ۱۶/۲۲ ۱٩‏ ۱۲۷ ۲۰/۶۴ ۲۸ ۳۷ شمالی 
واقع شده است (شکل ۱). ارتفاع متوسط این منطقه». ۲۶ متر پایین‌تر 
از سطح آب‌های آزاد است. میانگین سالیانه دما ۱/۷ سانتی گراد و 
میانگین بارندگی سالیانه آن ۱۳۵۶/۳۸ میلی‌متر است. میانگین حداکثر 
دما در گرم‌ترین ماه‌های سال (تیر و مرداد)" ۳۳ سانتی‌گراد و متوسط 
حداقل دما در سردترین ماه‌های سال (دی و بهمن)" ۰/۴ سانتی‌گراد 
بوده است ( ع بنا۱۱۵۲0 ۱۷۵۳۷۵ تلنمهگ ب‌مطمع‌هنطعد۲ 
6 ,۳۲۵۲101). 
داده‌های مورد استفاده 

در این تحقیق از تصاویر سال‌های ۲۰۰۰ و ۲۰۱۰ سنجنده‌ی 
ماهواره لندست ۵ و تصویر ۲۰۲۰ سنجنده‌ی آمآ) ماهواره 
لندست ۸ به منظور تهیه نقشه‌های کاربری اراضی استفاده شد 
(5.80۷عدنا10۲۵۲۰ص عمط تدع//:110). تصاویر, برمبنای میزان کمترین 
درصد ابرناکی و تقویم رشد محصول برنج در منطقه مورد مطالعه» از 
سایت زمین‌شناسی آمریکا دریافت شد (جدول ۱). همچنین, در این 
مطالعه از نقشه‌های کاربری اراضی تهیه شده در قالب شیپ فایلا 
که در دهه‌ی ۱۳۸۰ تهیه شده بود و نیز نقشه گوگل ارث ‏ برای تعیین 
نقاط کنترل زمینی و تصحیح هندسی تصاویر ماهواره‌ای استفاده 
گردید. به منظور پردازش تصاویر ماهواره‌ای» یافتن بهترین ترکیب 
باندی» تعیین نقاط تعلیمی و تهیه نقشه‌های کاربری مربوط به هر 
سال از نرم‌افزارهای 5.3 2۳۲۷۲ 3.3 11۷۲5 طاتقظ علعممن 
2 ۲ و 10.4 015 ۸۳ استفاده شد. 
پیش‌پردازش تصاویر 

در این مرحله, تصاویر از نظر هندسی» اتمسفری و رادیومتریک 
تصحیح شدند. ابتدا تصحیح هندسی بر روی تصاویر صورت گرفت. 
لازم به ذکر است تصاویر اخذ شده از سایت 7905 که با پسوند 
7 تهیه شده‌انده دارای تصحیحات هندسی می‌باشند. اما از آنجا که 
هدف این مطالعه مقایسه تصاویر چندزمانه می‌باشد برای انطباق 
کامل تصاویر سه سال تصویر سال ۲۰۲۰ به دلیل دقت بالاتر آن. 
به‌عنوان تصویر مرجع در نظر گرفته شد و دو تصویر ۲۰۰۰ و ۲۰۱۰ با 
روش تصویر به تصویر " با تصویر ۲۰۲۰ زمین‌مرجع شدند. همچنین, 
به دلیل اثرات جو بر روی میزان انرژی ثبت شده توسط سنجنده, 
برای تحلیل داده‌های تصویر می‌بایست اثر جو درصورت امکان تعدیل 
می‌شد. دراین تحقیق, به دلیل آن که ماهیت پردازش صورت گرفته» 
طبقه‌بندی بو اثرات جو بر روی نتایج تأثیر قابل‌توجهی نداشت. 
بنابراین از روش تصحیح اتمسفریک ساده تفریق پیکسل‌های تاریک" 


2- ۳0206 16 

طت2 000216 -3 
6 0۵ 19220 -4 
۵۵۵۵۵ 121[ -5 


مصری و همکاران 


استفاده شد. 


45000 


. پایش سطح زیر کشت برنج با استفاده از سیستم اطلاعات جغرافیایی ... ۵۱٩‏ 


000 


شکل ۱ - منطقه مورد مطالعه 
2 ناو عط 0۴ 100240100 ۲۱2۰1۰ 


جدول ۱- مشخصات تصاویر مورد استفاده در طبقه‌بندی 
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در حالت ایده‌ال» پدیده‌های تیره دارای تابنش صفر در همه طول 
موج‌ها هستند. در این روش فرض می‌شود که در هر باند از تصویر 
می‌توان پیکسل‌هایی مانند آب یافت که مقادیر بازتابندگی آن‌ها صفر 
یا نزدیک به صفر باشد. به این ترتیب اثر پخش جوی به صورت مقدار 
ثابت به پیکسل‌ها در هر باند اضافه می‌گردد. به همین جهت برای 
حذف خطاهای اتمسفریک باید مقدار ثابتی از ارزش پیکسل‌های هر 
باند کم شود (2012 ,۷2۵8۵262 > صدعدظ). 

پرای تشکیل تصویر با بیشترین قدرت تفکیک و نمایش بهتر 
تغییرات منطقه مورد مطالعه. از شاخص 01۳ برای تعیین بهترین 
ترکیب باندی استفاده شد. بیشترین میزان 011 برای تصاویر مربوط 
به هر سال, در جدول ۲ آورده شده است. 


۴۲ ۵ صباتام ۵ -1 


1۲ سنجنده 


52106 ماهواره 

5 لندسته۵ 
5 لندست۵ 
8 لندست۸ 


طبقه بندی تصاویر 

برای طبقه‌بندی تصاویر ماهواره‌ای» از طبقه‌بندی نظارت شده و 
به روش حداکثر احتمال استفاده شد. روش حداکثر احتمال» یکی از 
دقیق‌ترین روش‌های طبقه‌بندی در اکثر تحقیقات گزارش شده است 
۵0۵۲06 :1988 ملصهوعلانا ی مصهعامه مفصتا۳0۵) 
(1999 ,ول‌تهطت؟. در این روش, احتمال این که یک پیکسل بتواند 
به هر یک از کلاس‌های موجود تعلق یابد محاسبه می‌شود و سپس 
پیکسل به کلاسی که بیشترین احتمال را دارده اختصاص می‌یابد 
۱ ,1۳216027 ی ,۸20062۳02 رتتام‌ممهطاطاه ره 0۵02721 
212 برای محدوده‌ی مورد مطالعه هفت کاربری شامل اراضی 
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پهنه‌های ماسه‌ای 9 سایر اراضی تعیبن گردید. با استفاده از نرم‌افزار 
گوگل ارث و نقشه‌های موجود از کاربری اراضی منطقه مورد مطالعه با 


۰ ۲۰۲۰ تهیه شد. 


جدول - نتایج حاصل از شاخص 01۳ 
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مقدار شاخص ۷۵۱76 10067 
6229 
6293 
۱133-94 


اعتبارسنجی و ارزیابی صحت طبقه‌بندی 

جهت ارزیابی دقت طبقه‌بندی, از نمونه‌های تعلیمی‌ای استفاده 
شد که در روند طبقه‌بندی دخالت داده نشده‌بودند. برای هر بک از 
کلاس‌هاء ۷۵ درصد پیکسل‌ها به‌عنوان نمونه‌های آموزشی در 
طبقه‌بندی مورد استفاده قرار گرفتند و ۲۵ درصد مابقی پیکسل‌ها 
به‌عنوان نمونه‌های اعتبارسنجی در نظرگرفته شدند. ارزیابی دقت و 
صحت نقشه‌های طبقه‌بندی شده با استفاده از پارامترهای دقت کلی» 
ضریب کایاء دقت تولیدکننده , دقت کاربر» خطای اضافه و خطای 
حذف" انجام گرفت. دقت کلی» میانگینی از دقت طبقه‌بندی است و از 
نسبت تعداد پیکسل‌های درست طبقه‌بندی شده به کل پیکسل‌های 
طبقه‌بندی شده در تمام طبقات به‌دست می‌آید. یکی دیگر از 
پارامترهای دقت که از ماتریس خطا استخراج می‌شود ضریب کاپا 
است که دقت طبقه‌بندی را نسبت به یک طبقه‌بندی کاملا تصادفی 
محاسبه می‌کند (2016 ,150 6 ۷120067). مزیّت ضریب کاپا نسبت 
به دقت کلی» در استفاده از مقادیر حاشیه‌ای ماتریس خطا برای 
محاسبه دقت است. علاوه‌بر آن» دقت کی یک برآورد خوش‌بینانه 
استا و هبیشه دقت زا بالاتر از مقدار واقصی معصانسبه مس کنه. از 
طرفیء دقت کلی و ضریب کاپا با تمام طبقات سر و کار دارند و 
اطلاعاتی درباره هر کلاس نمی‌دهند. برای محاسبه دقت در هر 
کلاس مجزاء باید از پارامترهایی نظیر دفت کاربر و دفت تولیدکننده 
انتفاه قوس حول ماین مار خطای حافه اعافه کی بیان 
می‌گردد. خطای اضافه عبارت است از درصدی از پیکسل‌ها که متعلق 
به کلاس مورد نظر نبوده. ولی در آن کلاس قرارگرفته‌اند. خطای 
حذف» درصدی از پیکسل‌ها هستند که در اصل. متعلق به کلاس 
موردنظر بوده, ولی به اشتباه در کلاس دیگر قرار گرفته‌اند ( تداهحتلخر 
4 ,0 61). 
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ترکیب باندی 2705 سال ۱۷۵۵۲ 
1,4,5 2000 
1,4,5 2010 
4,5,7 2020 


ماتریس‌های خطای طبقه‌بندی در سال‌های ۰۲۰۰۰ ۲۰۱۰ و 
۰ شامل مقادیر دفت کلی» ضریب کاباء خطای اسافه خظای 
حذف» صحّت تولیدکننده و صحّت کاربر» به‌ترتیب در جداول ۲ تا۵ 
نشان داده شده‌اند. 

با قیخه رنه جداول 9۳۳ ۵ مشاهزه مشود که بیکا های طقه 
اراضی برنج از سایر طبقات به‌خوبی قابل تفکیک بوده و این امر ناشی 
از یکپارچگی, فشردگی و گستردگی کلاس اراضی بنج در منطقه 
است. به‌طور کلی هرچه که مزارع برنج متراکم‌تر بوده و فضای خالی 
در میان آن‌ها کمتر باشد. طبقه‌بندی تصویر از دقت بالاتری برخوردار 
خواهد بود. همچنین میزان دقت کلی برای سال‌های ۰۲۰۰۰ ۲۰۱۰ و 
۰ به‌ترتیب برابر با ٩۷/۵۹ ٩۸/۴۵‏ و ۹۸/۷۲ درصد است. میزان 
ضریبکاپا برای سال‌های ۰۲۰۰۰ ۲۰۱۰ و ۲۰۲۰ نیز به‌ترتیب برابر با 
۱ ۰/۹۷۰ و ۰/۹۸۴ است. نتایج حاصل از دقت کلی و ضریب 
کاپا در هر سال نشان از دقت بالای طبقه‌بندی دارد. بالاتر بودن 
ضریب کاپا در تصاویر جدیدتر به‌دلیل وجود نمونه‌های کنترل زمینی 
نزدیک‌تر از نظر زمانی به این سال‌ها و قدرت تفکیک بیشتر این 
قضاویا می‌باشه, ققه‌های کارت رای قیقر شک هاین ۲ 
تا ۴ به نمایش درآمده است. همچنین, مقادیر مربوط به مساحت هر 
یک از کاربری‌ها بر حسب هکتار و درصد در جدول ۶ ارائه شده است. 

در جدول ۷ میزان تغییرات هر یک از کاربری‌ها در طی سال‌های 
ما وا ۲۸۲۰ و جع قا +۲۲ مهرد پزرشی ق ار 
گرفته است. اراضی برنج در سال ۲۰۰۰ با مساحت ۱۱۰۸۰/۶۶ هکتار 
حدود ۵۸ درصد از کل کاربری‌های این منطقه را دربر می‌گیرد. 
مساحت این اراضی از سال ۲۰۰۰ تا ۲۰۱۰ کاهش ۴/۴۲ درصدی 
داشته است. همچنین در سال‌های ۲۰۱۰ تا ۲۰۲۰ با کاهش ۲/۶۴ 
درصدی رو به رو بوده است. در مجموع اراضی برنج از سال ۲۰۰۰ تا 
۰ با کاهش ۶/۹۴ درصدی همراه بوده و مساحت آن به 
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۹ هکتار رسیده است. افزایش تورم و قیمت زمین‌های 
زراعی» همچنین رشد جمعیت و ورود غیر بومیان به منطقه را می‌توان 
یکی از دلاثل تغییر کاربری اراضی زراعی به مسکونی دانست. 
ماخ رگا تیوه ی ال هی ما ۳ 
۸ درصد کاهش داشته است. اما بین سال‌های ۲۰۱۰ تا ۲۰۲۰ 
روند نزولی آن کاهش يافته و برابر با ۲۶/۲۶ درصد شده است. در 
مجموع بین دوره ۲۰۰۰ تا ۰۲۰۲۰ ۱۶۸۲/۲۵ هکتار معادل ۶۱/۳۲ 
درصد از اراضی جنگل‌های نیمه‌انبوه کاهش يافته است. این اتفاق 
به‌دلیل برش بی‌روبه درختان و تغفییر کاربری اراضی جنگل‌های 
نیم‌انوه به جنگل‌های تنک و مناطق مسکونی است. همچنین بین 
یال ای ۶ شالف جنگ های تک ۱/۳ درمند رگ 
داشته است که به دلیل کاهش وسعت جنگل‌های نیمه‌انبوه است. اما 
وا فا ها ایا اه ۷۱۲۶ در ها موه 
است که می‌تواند به دلیل تغییر کاربری به اراضی مسکونی باشد. در 


جنگل‌های تنک افزايش داشته است. بین سال‌های ۲۰۰۰ تا ۲۰۱۲۰ 
مساحت اراضی مسکونی ۲۲/۷۸ درصد رشد داشته است. همچنین در 
سال‌های ۲۰۱۰ تا ۲۰۲۰ مساحت آن با ۲۶/۴۸ درصد افزايش همراه 
بوده است. در مجموع مساحت مناطق مسکونی ۱۶۹۵/۸۴ هکتار 
معادل ۶۷/۹۴ درصد رشد داشته است. این رشد می‌تواند ناشی از 
افزايش جمعیت در منطقه» مهاجرت روستاییان به شهرها و در نتیجه 
تغییر کاربری اراضی زراعی و جنگلی به مناطق مسکونی باشد. 
پهنه‌های ماسه‌ای در سال‌های ۲۰۰۰ تا ۲۰۱۰ با کاهش مساحت به 
میزان ۵۸/۵۶ درصد رو به رو بوده است. همچنین بین سال‌های 
۰ تا ۰۲۰۲۰ ۷/۰۵ درصد کاهش داشته است. در مجموع مساحت 
کاربری پینه‌های ماسه‌ای بین سال‌های ۲۰۰۰ تا ۲۰۲۰ ۱۴۳/۱۵ 
هکتار معادل ۶۱/۴۸ درصد کاهش داشته است. این روند نزولی 


همراه با افرانقن ماخت وسازها در متاطق ساحلی همراه بوده انیت 


جدول ۳- ماتریس خطای طبقه‌بندی تصویر به روش حداکثر احتمال (سال ۲۰۰۰) 
(2000) ۵00 000حتام۱ن1 حصباصل‌جقهه عصتعیا معقصا ۵ تفه تمه ممتاهع1تعوهآه م1 -3 م1201 


پهنه‌های ‏ مناطقق ‏ مناطق 
کل سایر ماسه‌ای آبی مسکونی 
۵۱ 0067 90 ۵۵۲ مرن اند 
225 007 295 
108 0 0 0 0 
50 0 0 0 0 
55 0 0 0 0 
52 0 0 0 52 
28 0 0 28 0 
70 17 63 0 0 
89 89 0 0 0 
452 96 63 28 52 
0 10 0 0 
1.20 0 0 0 
9211 100 100 100 
100 90 100 100 
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55 50 108 کل 
ما1 
۵ ۳ خطای اضافه () 
(۶6) طمتوعتصصمت 
۳ 8 ۸ خطای حذف () 
(۶6) 100ووت0 
صحت تولید کننده () 
100 100 100 ۵۵ ۳۳۵0۱۵66۲ 
2 
صحت کاربر (/ 
100 100 100 بر (0) 
(0۶0 ۵۵۵۲۵0۷ ۲9687 
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جدول - ماتریس خطای طبقه‌بندی تصویر به روش حداکثر احتمال (سال ۲۰۱۰) 
(2010) 00مدظ 000ام۱ن1 حصباصل‌جقصه عصتعیا معقصصا ۵ تفص تمه 21100ع1تعوهآه 16 -4 101 
پهنه‌های مناطق مناطق جنکل‌های جنکل‌های اراضی 


کل سایر ماسه‌ای آبی مسکونی تنک نیمه انبوه برنج 
لها عطا0 5920 ۵۲ مانب 509 »عصعل-نصعد .. معز۳[ 
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(۵ 97.59) - ۸۵۵۷۵۲۵۸۷ 0۷۵۲۵11 دقت کلی (0.970) - 06111601601) 12002 ضریب کاپا 
بین سال‌های ۲۰۰۰ تا ۲۰۱۰ مساحت مناطق آبی ۳ درصد درصدی همراه بوده و در مجموع در اين بازه ۲۰ ساله مساحت مناطق 
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شکل ۲- نقشه کاربری اراضی سال ۲۰۰۰ 
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جدول ۵- ماتریس خطای طبقه‌بندی تصویر به روش حداکثر احتمال (سال ۲۰۲۰) 
(2020) ۵)800ظ 13۱۵11000 مصیاحی یمه عصتوا فعقحصا ۵۶ تمه تمه صمتاهع1/1دععاه م1 -5 1801 
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شکل 6- نقشه کاربری اراضی سال ۲۰۲۰ 
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جدول ۷- تغییرات مساحت هر یک از کاربری‌ها در سال‌های ۲۰۰۰ ۲۰۱۰ و ۲۰۲۰ 
0 220 2010 ,2000 1 1۷58 1240 مدع ۵۶ هعته2 عط ما مععصفط -7 13016 


2000-00 2010-00 2000-10 
(06) ۸۲69 (2ظ) وع۲ظ ۰۹9 (2ظ) وع۲ظ۸ ۰۹9 (2ظ) ۸۲62 
)6( )6( 
34.- 997 2.4 2-16 4.42 1( اراضی برچ 
1205 2:66 
جنگل‌های نیمه‌انبوه 
32- 65- 26.36 1 48 - 3- ععمعل- و۹ 
حا0169] 
1913 2918 7.26 8-- 29.03 9/0 ات 
٩02756 0655‏ 
694 19394 26048 9/8 22.78 91916 ۳ شییکوتي 
8 و« -]أتباظ 
4.1- 9.28- 25 29.20- 233 992 مناطق ابی 
07 ۱۷۷216۲ 
148 5 - 7.05 1- ۰-6 ۰ 136.34- پهنه‌های ماسه‌ای 
5 2۳00 
سا 
9611 6393 3.10- 42 99207 62,924 ی 


نتیجه گیری 

در تحقیق حاضر» به بررسی روند تغییرات کاربری‌های اراضی 
مختلف به‌ویژه سطح زیر کشت برنج در بخش کیاشهر شهرستان 
آستانه پرداخته شد. با توجه به نتایج به‌دست‌آمده میزان نرخ تغییرات 
کاربری برای کلاس‌های شالیزار و جنگل نیمه‌انبوه در بین سال‌های 
۰ تا ۲۰۲۰ کمتر از دوره ده ساله قبل از آن بود که این کاهش 
ذیربط مربوط دانست. اگرچه میزان کاهش این کاربری‌ها در دهه‌ی 
گذشته همچنان ادامه داشته است. در این میان, افزایش ۶۷/۹۴ 
نرخ تغییر کاربری در بخش مسکونی (با نادیده گرفتن کاربری سایر) 
نتیجه» لزوم توجه جدی به تغییرات کاربری در حوزه روستایی را بیش 
از پیش نشان می‌دهد. نکته دیگره رشد مساحت جنگل‌های تنک بین 
ال ها ها ها ۲ یی تن اب پی ال طای ۲۵ 


00134 


تا ۲۰۲۰ است که این نشان می‌دهد. در یک دوره ۱۰ ساله اقدام به 
تخریب جنگل‌ها شده و در دوره ده ساله بعد زمین‌های حاصل از این 
تخریب‌ها به کاربری مسکونی تغییر یافته است. در نتیجه لزوم توجه 
جدی ادارات و سازمان‌های منابع طبیعی به تخریب عمدی جنگل‌ها 
در طول زمان» با هدف تغییر کاربری آن‌ها به مسکونی و فروش آن 
ضروری به نظر می‌رسد. نتایج این مطالعه نشان می‌دهد که می‌توان 
از فناوری سنجش از دور و سیستم اطلاعات جفغرافیایی به‌عنوان 
روشی دقیق, قابل اطمینان و سریع برای مطالعه و بررسی» شناسایی و 
مقایسه کاربری‌های اراضی و تغییرات آن‌ها؛ به‌ خصوص کاربری 
زراعی و جنگل‌ها استفاده کرد. 


سپاسکزاری 

نویسندگان از حمایت مالی وزارت علوم» تحقیقات و فناوری و 
پارک علم و فناوری گیلان از اين پایان‌نامه در قالب کد اعتباری ۲۶- 
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